Search for: All records

Accretion discs properties should deviate from standard theory when magnetic pressure exceeds the thermal pressure. To quantify these deviations, we present a systematic study of the dynamical properties of magnetically arrested discs (MADs), the most magnetized type of accretion disc. Using an artificial cooling function to regulate the gas temperature, we study MADs of three different thermal thicknesses, hth/r = 0.3, 0.1, and 0.03. We find that the radial structure of the disc is never mostly supported by the magnetic field. In fact, thin MADs are very near Keplerian. However, as discs gets colder, they become more magnetized and the largest deviations from standard theory appear in our thinnest disc with hth/r = 0.03. In this case, the disc is much more extended vertically and much less dense than in standard theory because of vertical support from the turbulent magnetic pressure and wind-driven angular momentum transport that enhances the inflow speed. The thin disc also dissipates a lot of thermal energy outside of z/r = ±0.03 and a significant fraction of this dissipation happens in mildly relativistic winds. The enhanced dissipation in low-density regions could possibly feed coronae in X-ray binaries (XRBs) and active galactic nuclei (AGNs). Wind-driven accretion will also impact the dynamical evolution of accretion discs and could provide a mechanism to explain the rapid evolution of changing-look AGN and the secular evolution of XRBs. Finally, our MAD winds have terminal velocities and mass-loss rates in good agreement with the properties of ultrafast outflows observed in AGN.

 

Quasars show a remarkable degree of atomic emission-line broadening, an observational feature which, in conjunction with a radial distance estimate for this emission from the nucleus, is often used to infer the mass of the central supermassive black hole. The radius estimate depends on the structure and kinematics of this so-called broad-line region, which is often modeled as a set of discrete emitting clouds. Here, we test an alternative kinematic disk-wind model of optically thick line emission originating from a geometrically thin accretion disk under Keplerian rotation around a supermassive black hole. We use this model to calculate broad emission-line profiles and interferometric phases to compare to GRAVITY data and previously published cloud modeling results. While we show that such a model can provide a statistically satisfactory fit to GRAVITY data for quasar 3C 273, we disfavor it as it requires 3C 273 be observed at high inclination, which observations of the radio jet orientation do not support.

 

Magnetic reconnection is ubiquitous in astrophysical systems, and in many such systems the plasma suffers from significant cooling due to synchrotron radiation. We study relativistic magnetic reconnection in the presence of strong synchrotron cooling, where the ambient magnetization,σ, is high and the magnetic compactness,ℓ_B, of the system is of order unity. In this regime,e^±pair production from synchrotron photons is inevitable, and this process can regulate the magnetizationσsurrounding the current sheet. We investigate this self-regulation analytically and find a self-consistent steady state for a given magnetic compactness of the system and initial magnetization. This result helps estimate the self-consistent upstream magnetization in systems where plasma density is poorly constrained, and can be useful for a variety of astrophysical systems. As illustrative examples, we apply it to study the properties of reconnecting current sheets near the supermassive black hole of M87, as well as the equatorial current sheet outside the light cylinder of the Crab pulsar.

 

We explore reprocessing models for a sample of 17 hypervariable quasars, taken from the Sloan Digital Sky Survey Reverberation Mapping project, which all show coordinated optical luminosity hypervariability with amplitudes of factors ≳2 between 2014 and 2020. We develop and apply reprocessing models for quasar light curves in simple geometries that are likely to be representative of quasar inner environments. In addition to the commonly investigated thin-disk model, we include the thick-disk and hemisphere geometries. The thick-disk geometry could, for instance, represent a magnetically elevated disk, whereas the hemisphere model can be interpreted as a first-order approximation for any optically thick out-of-plane material caused by outflows/winds, warped/tilted disks, and so on. Of the 17 quasars in our sample, 11 are best-fitted by a hemisphere geometry, five are classified as thick disks, and both models fail for just one object. We highlight the successes and shortcomings of our thermal reprocessing models in case studies of four quasars that are representative of the sample. While reprocessing is unlikely to explain all of the variability that we observe in quasars, we present our classification scheme as a starting point for revealing the likely geometries of reprocessing for quasars in our sample and hypervariable quasars in general.

 

X-ray binaries exhibit a soft spectral state comprising thermal blackbody emission at 1 keV and a power-law tail above 10 keV. Empirical models fit the high-energy power-law tail to radiation from a non-thermal electron distribution, but the physical location of the non-thermal electrons and the reason for their power-law index and high-energy cut-off are still largely unknown. Here, we propose that the non-thermal electrons originate from within the black hole’s innermost stable circular orbit (the ‘plunging region’). Using an analytic model for the plunging region dynamics and electron distribution function properties from particle-in-cell simulations, we outline a steady-state model that can reproduce the observed spectral features. In particular, our model reproduces photon indices of Γ ≳ 2 and power-law luminosities of the order of a few per cent of the disc luminosity for strong magnetic fields, consistent with observations of the soft state. Because the emission originates so close to the black hole, we predict that the power-law luminosity should strongly depend on the system inclination angle and black hole spin. This model could be extended to the power-law tails observed above 400 keV in the hard state of X-ray binaries.

 

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of 10–100 for ≈60 min. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from particle-in-cell simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles, and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies γe ≲ 106 and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simplified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole, and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard (≲ 100 keV) X-ray energies.

 

Magnetically arrested accretion discs (MADs) around black holes (BHs) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here, we propose instead that it is mainly a dynamically important toroidal field – the result of dynamo action triggered by the significant but still relatively weak vertical field – that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotational instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (standard and normal evolution; strongly magnetized but not ‘arrested’) discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

 

Abstract We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87. 

ABSTRACT Infrared interferometry is a new frontier for precision ground-based observing, with new instrumentation achieving milliarcsecond (mas) spatial resolutions for faint sources, along with astrometry on the order of 10 microarcseconds (μas). This technique has already led to breakthroughs in the observations of the supermassive black hole at the Galactic centre and its orbiting stars, active galactic nucleus, and exo-planets, and can be employed for studying X-ray binaries (XRBs), microquasars in particular. Beyond constraining the orbital parameters of the system using the centroid wobble and spatially resolving jet discrete ejections on mas scales, we also propose a novel method to discern between the various components contributing to the infrared bands: accretion disc, jets, and companion star. We demonstrate that the GRAVITY instrument on the Very Large Telescope Interferometer should be able to detect a centroid shift in a number of sources, opening a new avenue of exploration for the myriad of transients expected to be discovered in the coming decade of radio all-sky surveys. We also present the first proof-of-concept GRAVITY observation of a low-mass XRB transient, MAXI J1820+070, to search for extended jets on mas scales. We place the tightest constraints yet via direct imaging on the size of the infrared emitting region of the compact jet in a hard state XRB.